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Abstract — In this paper, we consider quantum error correction 
over depolarizing cliannels witli non-binary low-density parity- 
check codes defined over Galois field of size 2^. The proposed 
quantum error correcting codes are based on the binary quasi- 
cyclic CSS (Calderbank, Shor and Steane) codes. The resulting 
quantum codes outperform the best known quantum codes and 
surpass the performance limit of the bounded distance decoder. 
By increasing the size of the underlying Galois field, i.e., 2'', the 
error floors are considerably improved. 

Index Tenns — LDPC code, non-binary LDPC codes, belief 
propagation, Galois field, iterative decoding, CSS codes, quantum 
error-correcting codes 

I. Introduction 

In 1963, Gallager invented low-density parity-check 
(LDPC) codes fl]. Due to the sparseness of the code represen- 
tation, LDPC codes are efficiently decoded by the sum-product 
algorithm. By a powerful optimization method density evo- 
lution [2 1, developed by Richardson and Urbanke, messages 
of sum-product decoding can be statistically evaluated. The 
optimized LDPC codes exhibit error correcting performance 
very close to the Shannon limit f3l. Recently, LDPC codes 
have been generalized from a point of view of Galois fields, 
i.e. non-binary LDPC codes are proposed. Non-binary LDPC 
codes were invented by Gallager [l]. Davey and Mac Kay Q 
found that non-binary LDPC codes can outperform binary 
ones. 

Quantum LDPC codes, which are quantum error-correcting 
codes, have been developed in a similar manner to (classical) 
LDPC codes. By the discovery of CSS (Calderbank, Shor and 
Steane) codes [5|, |6| and stabilizer codes [7|, the notion 
of parity-check measurement, which is a generalized notion 
of parity-check matrix, is introduced to quantum information 
theory. In particular, a parity-check measurement for a CSS 
code is characterized by a pair of parity-check matrices which 
satisfy the following condition: the product of one of the pair 
and the transposed other is subjected to be a zero-matrix. 

Quantum LDPC codes were first introduced by Postol in |8). 
The above CSS constraint on the parity-check matrices makes 
the design of the quantum LDPC codes difficult. MacKay et 
al. proposed the bicycle codes |9| and Cayley graph based CSS 
codes ifTOl . In ifTTl . Poulin et al. proposed serial turbo codes 
for the quantum error correction. These codes can be decoded 
by an efficient iterative decoder. To the best of the authors' 
knowledge, these codes |9l, IfTOl . IfTTl are the best known 
quantum error correcting codes among efficiently decodable 
quantum LDPC codes so far. In |12|, Hagiwara and Imai 
proposed a construction method of CSS code pair that has 



quasi-cyclic (QC) parity-check matrices with arbitrary regular 
even row weight L > A and column weight J such that 
i/2 > J > 2. However, the resulting codes do not outperform 
the codes proposed by MacKay el al. |9|, lITOl . 

Generally, LDPC CSS codes tend to have poor minimum 
distance. The minimum distance of an LDPC CSS code is 
upper-bounded by the row weight of the parity-check matrix. 
This is due to the dual and sparse constraint on the parity- 
check matrices. When the LDPC CSS codes are used with 
large code length, the poor minimum distance leads to high 
error floors. Therefore, it is desired to establish the construc- 
tion method of quantum LDPC codes with large minimum 
distance. We should note that it is important to study quantum 
LDPC codes with large minimum distance which grows with 
code length fT3l for constructing quantum LDPC codes with 
vanishing decoding error probability. 

Non-binary LDPC codes are defined as codes over GF(2p) 
with p > 2. The parity-check matrices of non-binary LDPC 
codes are given as sparse matrices over GF(2^). In this paper, 
we investigate non-binary LDPC codes for quantum error 
correction. It is empirically known that the best classical non- 
binary LDPC codes have column weight J — 2 from a point 
of view of error-correcting performance |[T4l . Moreover, due to 
the sparse representation of non-binary parity-check matrices 
of column weight J = 2, the non-binary LDPC codes are 
efficiently decoded by FFT-based sum-product algorithm |15|. 

In this paper, we propose a construction method of a 
binary CSS code which can be viewed also as a pair of non- 
binary LDPC codes. More precisely, the proposed construction 
method produces a binary code pair (C, D) such that C D D^, 
and C and D are also defined by non-binary sparse parity- 
check matrices over GF(2'°) of column weight J = 2. 
This satisfies the constraint of CSS codes. To this end, we 
first construct PJ x PL binary QC parity-check matrix pair 
{Hc,Hd) with column weight J = 2 and row weight L 
such that HcHJ) = by the method developed in lfT2l . 
Solving some linear equations on Z2P-1, we get PJ x PL 
non-binary parity-check matrix pair (-ffr,^A) with column 
weight J = 2 and row weight L such that HrH\ = 0. It is 
known that a natural Unear map from GF(2p) to GF(2)p^p is 
given so that through this map, the non-binary LDPC matrix 
pair {HytH/\) can be viewed as a binary LDPC matrix pair 
[HctHd) such that HcHj^ — 0. Numerical experiments 
show the resulting CSS codes outperform the best known 
quantum error correcting codes and surpass the performance 
limit of the bounded distance decoder. By increasing the size 
of the underlying Galois field, i.e., 2^, the error floors are 
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Fig. 1. An example of binary {J = 2,L = 6, P = 7)-QC parity-check matrix pair {He, Hu) constructed by the method in Theorem [T] with a = 2 and 
T = 3. It holds that HcHj^ = 0. For any row m' of //j5 



considerably improved. 

The rest of this paper is organized as follows. Section |ll] 
describes the construction method of a non-binary twisted 
LDPC parity-check matrix pair {HG,Ha) of column weight 
J = 2. Section |III] describes the decoding algorithm of the 
binary twisted code pair {C,D). Section HV] demonstrates the 
decoding performance of the proposed codes. 



II. Construction of Non-Binary Matrix Pair with 
Column Weight 2 

In this section, we will construct a binary code pair (C, D) 
defined by orthogonal parity-check matrices He and Hd, 
where we call two matrices X and Y orthogonal if XY~^ = 0. 
Let pN qubits be the code length. The binary codes C and D 
are designed in such a way that C and D can be represented 
also by non-binary sparse parity-check matrices over GF(2p) 
of column size N and column weight J = 2. To this end, we 
start with binary QC matrices and extend them to matrices 
over GF(2'°). The following is the outline of construction 
procedure. 

A) Choose integers J = 2,P and L such that PL = N. 
Construct a pair of orthogonal PJ x PL quasi-cyclic 
binary matrices He and Hd by following the procedure 

B) A pair of PJ x PL matrices Hr and Ha over GF(2p) are 
constructed by replacing non-zero entries of the matrices 
He and Hjj with elements of GF(2P). The orthogonality 
condition HpHj^ — imposes a set of linear constraints 
on the non-zero entries of Hp and Ha, that can be solved 
by Gaussian elimination. 

C) The entries of Hr and Ha are mapped to elements 
of GF(2)P^P that preserves addition and multiplication. 
These properties of the mapping imply that the resulting 
matrices He and Hd remain orthogonal. 

The proposed code can be viewed as a concatenated code with 
an inner LDPC code and an outer code of rate 1 and length 
p. Indeed, the code obtained at step A) is a quantum LDPC 
code |16|. Steps B) and C) have the effect of replacing each 
qubit by p qubits while preserving the orthogonahty property. 
Each step in the procedure is explained in the following 
sub-sections. 



A. Quasi-Cyclic Binary Construction 

Let He and Hd be PJ x PL binary parity-check matrices 
defined as follows: 

He '■= iIicj/))o<3<j,o<i<L, 
Hd '■= iIidjj))o<j<j,Q<e<L: 
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We refer to such matrices as (J, L, P)-QC matrices. 

Hagiwara and Imai proposed |12| the following method for 
constructing a QC parity-check matrix pair [HeiHo)- In the 
original paper [121, the construction method is more flexible 
about the row size of the matrices, i.e.. He and Hd can have 
different row sizes. For simplicity, in this paper, we focus on 
He and Hj:) with the same row size JP. 

Theorem 1 (|16, Theorem 5.2]). Define Z*p ■- {z e Zp \ 

3a G 1jp,za = 1}, and ord(cr) :— niin{?77, > | cr™ — 1}. 
For integers P > 2, J, L,0 < a < P and < t < P such 
that 



a,T e Zp, 

L/2 = ord(cr), 

1 < J < ord(o-), 

ord(a) ^ #Z*p, 

1 - 0-' £ Zp for all 1 < j < ord(cr), 

r ^ l,cr,cr^,...,cr 



ord(CT)-l 



(1) 

(2) 



(3) 
(4) 



let He and Hd be two (J, L, P)-QC binary matrices such 
that 



He = {l{cjj))o<j<j.i3<i<L, 

Hd = {I{d-j,i))o<j<j,o<i<L, 

„-3+e 0<i<L/2 
L/2<1< L, 






-j+e 



0<i<L/2 
L/2<£ < L 
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(6) 
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Fig. 2. An example of non-binaiy matrices Hr = {-fm,n)o<m<M,o<n<N and Ha = {Sm.,n)o<m<M,o<n<N over GF(2*) such that HrH\ = 
with A/ = 14 and N = 42. Each non-zero entry is represented as the hexadecimal number of log^{'ym,n), where o is a primitive element such that 
a* + a + 1 = 0. E.g., a" and a^^ are represented as and b, respectively. 



then it holds that HcHJ) — and there are no cycles of size 
4 in the Tanner graph of He and Hd. 

From Theorem [T] we obtain two JP x LP binary matrices 
He and Ho such that HcHj, ~ and the Tanner graphs of 
He and Hd are free of cycles of size 4. We give an example. 

Example 2. With parameters J — 2,L — 6,P~7,<7^2 
and T = 3, from Theorem |7] we are given a JP x LP binary 
matrix pair (He, Hd) such that Hc'Hj) = as follows. 

^ _fl{l) m /(4) /(3) /(6) /(5)^ 
''^-l,/(4) /(I) /(2) /(5) /(3) /(6), 

fr _fm m /(I) /(6) /(3) /(5)^ 

""""-[lil) /(4) /(2) /(5) /(6) /(3), 

The binary representation of these matrices are given in Fig. [7] 
It can be verified that there are no cycles of length 4 in Tanner 
graphs of He and Ho. 

We observe a fundamental property of He and Ho as 
follows. 

Observation 3. Let m' :— 5. The m'-th row of Ho has non- 
zero entries at the uq = 2, 712 = 7, n^ = 20, ni ~ 25, n^ — 29 
and ns — 38-th columns. Denote the set of these indices by 
N{m') := {no, ... ^n^}. Note that the index starts from 0. 
These non-zero entries are represented in thick font in Fig. [7] 
At each of these 6 columns in N{m'), He has J = 2 non-zero 
entries at 

( 1,2), ( 5,7), (2,20), ( 1,25), ( 2,29), ( 5,38), 
(12,2), (13,7), (11,20), (13,25), (12,29), (11,38). 

Let these positions be denoted, respectively, by 

{mo, no), (m2,n2), (to4,'^4), (mo,ni), (7714,^5), (nT'2,"3), 
(m5,no), (mi,n2), (7773,714), (7774,771), (7775,775), (nis^n^). 

We denote the set of these positions by E{m'). In the Tanner 
graph of He, those non-zero entries from (ttiq, 779) to (7775, no) 
consist of a cycle of size 2L = 12. Along this cycle, the row 
index moves back and forth between the above and below 
of He. Indeed, < 7772i < P < m2i+i < 2P for i = 
0, . . . , L/2 — 1. On the other hand, the column index moves 
back and forth between the left and right of He. Indeed, < 
772i < LP/2 < 772i+i < LP for i ^ Q, . . . , L/2 - 1. 

We claim that this observation is general for any 777'-th row 
of De- The following lemma will be a key ingredient for 
constructing non-binary matrices i^A and Hy in Section ITl-BI 



Lemma 4. Let He , Ho be the two (2, L, P)-QC binary 
matrices dealt in Theorem\l\ Let N{ni') := {no, . . . ,77^-1} 
be the support of the m'-th row of Ho. To be precise, 

N{m') = {no, . . . ,77i_i} = {0 < 77 < LP I d.m',n ^ 0}. 

Let Eim!) be the set of non-zero entry positions in He whose 
column is in N{m'). To be precise. 



E{m') := {(777,77) I c„^„ ^ 0,77 e iV(?77')}. 



(7) 



In this setting, in the Tanner graph of He, for any m' — 
0, . . . , i — 1, the L variable nodes corresponding to the column 
index in N{m') and the L adjacent check nodes form a cycle of 
length 2L. In other words, there exist L distinct mo, ■ • ■ , ni^^i 
and L distinct no, ... , 77i_i, such that 



{im2i-l,n2i), (77721, 772i), (»772i, 772i+l) , (7772i+l , 772i+l) 
I < i < L/2} = £:(777'), 

where we denote m^i := m^^i and 777^+1 := 777i. 



(8) 



Proof. Since J = 2, it follows that He and Ho comprise 
2x L sub-matrices of size P x P. For simplicity, we focus on 
the 777'-th row chosen from the upper half rows of Ho, ie., 
< 777' < P. The proof for P < m' < 2P is essentially the 
same. 

First, we will clarify the support N{m') of the 777'-th row 
of Ho- From (|6]), the upper half of Ho can be written by L 
sub-matrices as follows. 



■/(- 



-(L/2-1) 



)l|/(- 



-L/2\ 



■H- 



-(L-1) 



)) 



For each sub-matrix, the 777'-th row has the only one non- 
zero entry. For < i < P/2, let 772i be the column index 
of such a non-zero entry in the i-th sub-matrix. Similarly, 
let 772i+i be the column index of such a non-zero entry 
in the ([— 7]l/2 + i/2)-th sub-matrix. Obviously, N{m') = 
{no, . . . , 77,L-i} and 770, . . . , 77l_i are distinct. Note that the 
773-th, . . ., 77L-3-th and 77i-th columns are in the i/2 right 
half (i - l)-th, . . . , (i/2 + l)-th and (P/2)-th sub-matrices 
of Ho, respectively. The 777'-th row of I{x) has the only non- 
zero entry at the {x + 777')-th column. From this, it follows 
that we can rewrite 



"2j = [-ra ' +m']p + iP, 

n2^+l = [-(jl-'l^/^ + 7T7']p + i[-l]L/2 



L/2)P, 



(9) 
(10) 
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Fig. 3. An example of binaiy pM X pAf matrices Hq and Hu such tliat HqHJ^ = witli p = 4, M = 14 and TV = 42. Non-zero entries are represented 
in black. The codes have many cycles of size 4 as binary codes. On the other hand, the codes have no cycles of size 4 as non-binary codes. 



where we define [x]t G Z for x G Z as < [x]t < t such that 
[x]t — X (mod t). 

Secondly, we will prove ([8]l by clarifying the structure of 
E{m'). Define 



TO2i 






■^'' + m']p 



P. 



We claim that {m2i-i,n2i), {m2i,n2i) £ Eim') for < i < 
L/2. Consider the non-zero entry in the (j, i)-th sub-matrix 
for j = 0, 1, and < i < L/2. Let the position be denoted by 
(to, n2i)- From (|5]l and i < L/2, it follows that the 77i-th row 
of He has non-zero entries at the ([cr~^+* + m]p + iP)-th 
column. Therefore, it follows from (|9]l that (to, n2i) £ E(rn') 
if and only if 



,-j+« 



m\p 



]p + iP= [-Tu- 



rn \p 



iP, 



for J = 0, 1. Thus, we conclude that to ~ m2i-i and to = m2i 
for j = 0, 1, respectively, which proves the claim. In a similar 
manner, it can be shown that {'m2i,n2i+i), {m2i+i,n2i+i) G 
E{m') for < i < L/2, where we denoted to^+i := toi. 
These prove (HJ. 

Finally, we show that the 2L non-zero entries in E{m!) 
consist of a cycle of length 2L. To this end, it is suffice to show 
that TO-o, TO,2, . . . , 771L-2 are distinct, and ?7ii, TO3, . . . , tti^-i 
are distinct. Assume TO2i ~ m2i' for i ^ i'. Then we have 0-' + 
TCT^' + to' = (t* + Tcr^' + to' (mod P). Some calculations 
reveal that (1 - cr*'^')(cr* - rcr~*') = (mod P). From Q 
it follows that ct* — rcr^' = (mod P), which contradicts 
with dU. Hence, we conclude i — i'. In the same manner, we 
can show toi, TO3, . . . , ttil^i are distinct. D 

B. Non-binary matrix pair construction 

Define M := JP and N := LP. So far, we obtain 
orthogonal matrices He = {cm,n)o<m<ALQ<n<N and Hd = 
idm,n)o<m<M,o<n<N whose Tanner graphs are free of cycles 
of size 4. In this sub-section, we will construct orthogonal 
non-binai-y M x N matrices iJp = (7m,n)o<m<M,o<n<JV and 
Ha = {Sm,n)o<ni<M,o<n<N over GF(2P) such that 7„j,„ 7^ 
iff Cmji ^ and 6m.7i 7^ iff (i„j „ ^ 0. Obviously, the 
Tanner graphs of Hr and Ha are free of cycles of size 4. We 
will determine the non-zero entries of Hp and Ha such that 
HtHJ^ = in the rest of this sub-section. 



For HrHj^ = 0, it is required that the TO,'-th row of Ha is 
in the null-space of Hr for each < to' < JP. From Lemma 
m this is equivalent to 



"fmo,no lmo,n-i 



fmt-i,no 



1mL-2,nL-2 lmL-2,nL-l 
Tmi,-i,ni-i 



^m' ,nQ 



(11) 

for each < to' < JP. In order to find the non-zero entries of 
Hr and Ha, this equation needs to have non-trivial solutions, 
i.e., the determinant of the left matrix, denoted by r,„', in (fTTT i 
is 0: 



det(r™/) = 7„ 



•7- 



mL-i.riL-i 



(12) 
Divide E{m') in ^ into two parts as in the proof of Lemma 

El 

E{m') =Ei{m')UE2{m'), 

Ei{m') := {(TOo,no),(TOi,ni),...,(TOL_i,nL-i)}, 
E2{m) := {(TOo,ni),(TOi,ri2),...,(TOL_i,no)}. 
Then ( fT2] i can be transformed to 

n ^™-" n 7™'„=i- (13) 

{m,n)£Ei{m') {m,n)£E2{m') 

For a^ G GF(2p), define log„(a=") := a; (mod 2? - 1). Then 
logo, is well-defined. The equation above is equivalent to the 
following linear equation over Z2P_i. 

Yl l0ga7m,n- Yl loga 7m,n = 0. (14) 

(rn,n)G-Ei(m') {7n,n)^E2{m') 

Thus, we have JP linear equations over Z2P-1 for to,' = 
0, . . . , JP ~ 1. Solving these linear equations by the Gaussian 
elimination, we get the candidate solution space of the non- 
zero entries of Hp such that (fT4l i holds for to' = 0, . . . , JP—1. 
Picking non-zero entries of Hp randomly from the candidate 
solution space and solving ( fTTT i. we obtain non-zero entries of 
Ha. We give an example. 

Example 5. Using He and Hd given in Example^ we get an 
M X N non-binary matrix pair (Hp, Ha) over GF(2P) such 



that HrHl = with M = JP ^ 14 and N ^ LP = 42. entry of HrHj^. Then, for any < m < M and < n < TV, 
The resulting {Hr, H^) is depicted in Fig. |2] 



This construction can be viewed as picking (HtjHa) 
randomly from {{Hr, Ha) \ HtH\ = 0}, where Hr and Ha 
are constrained to have non-zero entries at the same positions 
as He and He,, respectively. Since He and He, is equivalent 
with some column permutation [12], the construction has 
symmetry for Hy and Ha- This symmetry leads to almost 
the same decoding performance which will be observed by 
computer experiments in Section HV] 



C. Binary Quasi-Cyclic CSS LDPC Codes 

So far, we obtained M x N sparse non-binary GF(2p) 
parity-check matrices Ha and Hy, where N := PL and 
M :— PJ. It is known that non-binary codes have the binary 
representation of their parity-check matrices. In this section, 
we show that two parity-check matrices Hy and Ha over 
GF(2'') such that HyHJ^ = can be represented by two 
binary matrices He and Hd such that HcHj) — 0. 

Let GF(2'') has a primitive element a with its primitive 



polynomial Tr{x) = J2i=o "" 



cP. It is known fl7| that 



the following map A from GF(2p) to GF(2)p^p is bijective 
and its image is isomorphic to GF(2'') as a field by sum and 
multiple as matrices. 

GF(2P) 3 a' ^ A{a') := A{ay G GF(2)p^p, 
A(0) = 0, 

"0000 TTo ■ 

10 TTl 

Aia):= 10 n, 

1 7rp_i 
Moreover, it holds that 

A{a')via^) = via'+'), 

where a' = Ej^o Oj"^ ^ GF(2p), 

and v{a') := (ao, . . . , flm-i)^ E GF(2)? 



Furthermore, with an abuse of notation we define A{v{a^)) := 

Fact 6. Let Hy and Ha be matrices over GF(2p)*^^^ and 
let He and Hu be two matrices over GF(2)p*^^p^ such that 

Hy = (7m,n)o<m<Af,0<n<Af7 
Ha = {Sm,n)o<m<M,0<n<N , 
He = (^(7m,n))o<m<M,0<n<JV, 
Hd = {A {S„i^n))o<m<M,0<n<N ■ 

Then, it holds that if HyH\ = 0, then HeHjj = 0. 

Proof: Let {HeHj))^ n be the (m, n)-th p x p binary 
sub-matrix of HeHj^, and let (-ffr^X)™." ^e the (TO,n)-th 



{HeHj))„i,n = J2k=Q Mlrn,k)A{Sk,n) 
= Y.k=0 Mlrn.kSn.k) 
= MJ2k=0 Im.kSn.k) 

= A{{HyHI\^,-,) = A{Q) = 0. 



Example 7. Using Hy and Ha given in Example\5\ we get a 

pMxpN binary matrix pair {He, Hd) such that HeHj, — 
with p = 4:,pM ^ pJP = 56 and pN = pLP = 168. The 
resulting {He, Hjj) is depicted in Fig. \3\ 

III. Decoding Algorithm 

In this section, we describe the decoding algorithm for the 
CSS code pair {C, D) constructed by the proposed method 
in Section |II] and III-CI The decoding algorithm is based on 
the decoding algorithm of classical non-binary LDPC codes 
[ISil. The input of the decoding algorithm is the syndrome. 
We assume the depolarizing channels |9 Section V] with 
depolarizing probability 2/m/3, where /,„ can be viewed as 
the marginal probability for X and Z errors. 

Let Af X iV be the size of the non-binary parity-check matrix 
Hy over GF{2p). The code length is pN qubits. We deal 
with a p-bit sequence as a non-binary symbol which is simply 
referred to as symbol. Moreover, we deal with the symbol 
interchangeably as a symbol in GF(2p). 

Note that the channel is the normal depolarizing channel. 
We assume the decoder knows the depolarizing probability 
3/m/2. For each row m = 1, . . . , Af in Hy, let Nm be the set 
of the non-zero entry indices in the m-th row. To be precise, 
Njn '■— {n I 7m. n 7^ 0}. The decoder is given the syndrome 
symbols s^ G GF(2)^ for m = 1, . . . ,M. To be precise, the 
decoder does not know the flipped qubits but their syndromes: 



Y^ A{'J,n,n)y^ 



(15) 



neN„ 



where A is the isomorphism defined in Section III-CI and y G 
GF(2)^ is a p-bit sequence corresponding to the n-th p-qubit 
sequence of flipped pN qubits. 

For simplicity, we concentrate on the decoding algorithm for 
C, since the decoding algorithm for D is given by replacing 
F with A, and A{-) with A'^ {■) in the following algorithm. 
The decoding algorithm of C 
initialization : 

For each column n = 1, . . . ,N in Hy, let A'/„ be the set of 
the non-zero entry indices in the n-th column. To be precise, 
Af„ := {m I 7m, n y^ 0}. For each column n in Hy for n = 
1, . . . ,N, calculate the initial probability p„ (e) as follows. 

p'^He) = Pr(e„ ^ e|y„ = 0) = /,^(^'(1 - /m)^"^"'-' 

for e G GF(2)P, where /m is the flip probability of the 
channel and Wfi{e) is the Hamming weight of e. For each 
column n = l,...,iV in Hy, copy the initial message 
Pnm — Pn E [0, 1]^'' for TTi G M„. Set the iteration round as 
£:^0. 



horizontal step 



M 



Each row m has L incoming messages pnm for v G N„i. 
The m-th row sends the following message qinn G [0,1]^'' 
to each column n G Nm- 



Pli(e) = pl!liAh-l)e) for e G GF(2)^ 

ymn — -"-s^ V^ 



^ n'm ' 



n'eW„^\{n} 

^n^He) - ?trH^(7«™k) for e G GF(2)^ 



(16) 



(17) 



where Is is a probability on GF(2)p such that Is (e) — 1 
for e = s„ and otherwise, and q\ ® q2 G [0, l]^*" is a 
convolution of qi G [0, 1]^ and 52 G [0, 1]^ .To be precise. 



{qi®(l2){e) 



E 9i(/)92(5)foreGGF(2)P. 



/,£eGF(2) 



The convolutions are efficiently calculated via FFT and IFFT 
ifTSl . |fT9l . Increment the iteration round as ^ := ^ + 1. 

vertical step : 

Each column n = 1,...,N in TJa has J = 2 non-zero 
entries. Let A/„ be the set of the column indices of the 
non-zero entry. The message piim G [0,1]^ sent from n to 
m G M„ is given by 

m'eAf„\{m} 

where ^ is the normalization factor so that 

L^e^GF(2)'pPn.m{S.) = I- 

tentative decision : 

For each n = 1,...,N, the tentatively estimated u-th 

transmitted symbol is given as 

el^) = argmaxpl°)(e) TT g^L(e). 

If (eg, . . . , e^) has the same syndrome as {s_i, . . . , s^/) which 
is defined in ( fTSb . in other words, for m = 1, . . . , iP, 



n£Nrr. 



^(7mn)ii'^=l™eGF(2)P 



for all c = 1, . . . , M, the decoder outputs (cq, . . . , e^) as the 
estimated error. Otherwise, repeat the latter 3 decoding steps. 
If the iteration round i reaches a pre-determined number, the 
decoder outputs FAIL. 

Note that, in this algorithm, the correlations between X 
errors and Z errors are neglected. In [9^ Section VI, C] MacKay 
et al. used the knowledge about the channel properties for 
decoding, which improved the decoding performance. The 
most complex part of the decoding is the horizontal step, 
which requires 0{Nq\og{q)) multiplications and additions 
when calculated via FFT, where q = 2^. 

IV. Numerical Result 

In this section, we demonstrate the proposed CSS code pair 
decoded by the algorithm described in the previous section. 
The proposed CSS code pair (C, D) is constructed as follows. 
First, by Theorem [T] construct JP x LP binary matrices He 



and Hf) with ,/ = 2,L, and P. Secondly, by the scheme 
described in Section [III construct JP x LPw non-binary 
matrices Hr and H^ over GF(2p). Finally, by the scheme 
described in Section lTl-CI we have pJPxpLP binary matrices 
He and Hd. Thus, we obtain C and D are defined by the 
parity-check matrices He and Ho, respectively. Note that C 
and D can not only be viewed as binary codes defined by He 
and Hd but also be viewed as non-binary codes defined by 
Hr and Hj\. The code length of the proposed CSS code is 
given as n = pLP qubits or equivalently LP symbols. The 
quantum rate Rq of the proposed CSS code is given as 

Rq^I- 2J/L. 

Fig. |4] shows the block error probability of the constituent 
codes C and D of the proposed CSS code pair (C, D) over 
the depolarizing channel with marginal flip probability f^ of 
X and Z errors. Parameter are chosen J = 2, L = 6,8 and 
14 for Rq — 1/3, 1/2 and 5/7, respectively. The depolarizing 
probability is given by 3/„i/2. The correlations between X 
errors and Z errors are neglected. Due to the symmetry of 
construction of C and D, the block error probability of the 
constituent codes C and D are almost the same, hence we plot 
the block error probability of either C or D. It is observed 
that for fixed q = 2^ and i?Q, the codes with larger code 
length tend to have higher error floors. This is due to the fact 
that the proposed codes have poor minimum distance which 
is upper-bounded by pL. The error floors can be improved 
by using larger p, i.e., larger field GF(2p), which leads to 
the requirements of more complex decoding computations 
0[Nq\og{q)), where q = 2^. 

Fig. |5] compares the proposed quantum codes with the 
best quantum codes so far The horizontal axis is the flip 
probability at which the block error probability of one of the 
constituent classical code is 0.5 x lO^"'. The vertical axis is 
the quantum rate i?Q of quantum codes. Since the proposed 
CSS codes have constituent classical codes C and D of the 
same classical rate Rq ~ 1 — J/ L, the quantum rate i?Q is 
given as i?Q — 2Rc — 1 = 1 — 2J/L. It can be seen that the 
proposed codes outperform the state-of-the-art codes. In fact, 
the proposed codes surpass the BDD curve which is the limit 
of the bounded distance decoder, while the other codes fall 
inside the BDD curve. 

V. Conclusion 

We proposed a novel construction method of CSS codes. 
The resulting CSS codes can be viewed as non-binary LDPC 
codes over GF(2p). Due to the sparse representation of 
the parity-check matrices, the proposed codes are efficiently 
decoded. The simulation results over the depolarizing channels 
show that the proposed codes outperform the other quantum 
error correcting codes which exhibited the best decoding 
performance so far The error floors are lowered by increasing 
the size of the underlying Galois field, i.e., 2^. 
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Fig. 5. The performance of the proposed CSS code pair (C, D) compared with the best CSS codes so far from f9], flOl and fill over the depolarizing channel 
with marginal flip probability /m of X and Z errors. Each point is plotted at which the block probabilities of both two constituent codes are 5 X 10~^. The block 
probability of the entire CSS code is 1 — (1 — 5 X 10^^)^ ~ 10~^. The Shannon Hmit of the depolarizing channel: _Rq = 1 — ?i(3/m/2) — 3/m/2 log2(3), 
where h{-) is the binary entropy function. The curve labelled S2 is the achievable quantum rate if the correlations between X errors and Z errors are neglected: 
Rq = 1 — 2/i(/m). The curve labeled BDD is the performance limit when the bounded distance decoder is employed and the correlations between X en'ors 
and Z errors are neglected. _Rq = 1 — 2/i(2/m)- The code length is n qubits. The proposed codes are defined over GF(q). 



